* Reproduces results in Lake, HIR, Tables A.2 and 5.3* use "Lake_HIR_ICB.dta" * Updated 8.6.18* Descriptive statistics for ICB data, Table A.2sum us_SH1995 us_indally us_indexmilperpop1995 us_EH1995 us_indexexchangerate1995 us_indexreltradep1995 duration geo_salience violence usinv* MULTINOMIAL LOGIT TEST*drop b*		*mlogit usinv us_SH1995 us_EH1995 duration violence geo_salience if usactor~=2, base(1) robust		*estimates store a1		estsimp mlogit usinv us_SH1995 us_EH1995 duration violence geo_salience if after_week==1, base(1) robust		*estimates store a2		setx meansetx us_SH1995 0 simqi, prval(2) genpr(start)setx us_SH1995 75simqi, prval(2) genpr(end)gen diff=end-startsumqi diff, level(95)drop start end diff		 setx meansetx us_SH1995 0 simqi, prval(3) genpr(start)setx us_SH1995 75simqi, prval(3) genpr(end)gen diff=end-startsumqi diff, level(95)drop start end diffsetx meansetx us_SH1995 minsimqi, prval(2) genpr(pr1)setx us_SH1995 p75simqi, prval(2) genpr(pr2)gen relrisk = pr2/pr1sumqi relriskdrop pr1 pr2 relrisksetx meansetx us_SH1995 minsimqi, prval(3) genpr(pr1)setx us_SH1995 p75simqi, prval(3) genpr(pr2)gen relrisk = pr2/pr1sumqi relriskdrop pr1 pr2 relriskdrop b*		estsimp mlogit usinv us_indally us_indexmilperpop1995 us_indexexchangerate1995 us_indexreltradep1995 duration violence geo_salience if after_week==1, base(1) robust		*estimates store a1		 setx meansetx us_indally 0 simqi, prval(2) genpr(start)setx us_indally 75simqi, prval(2) genpr(end)gen diff=end-startsumqi diff, level(95)drop start end diff		 setx meansetx us_indally 0 simqi, prval(3) genpr(start)setx us_indally 75simqi, prval(3) genpr(end)gen diff=end-startsumqi diff, level(95)drop start end diffsetx meansetx us_indally minsimqi, prval(2) genpr(pr1)setx us_indally p75simqi, prval(2) genpr(pr2)gen relrisk = pr2/pr1sumqi relriskdrop pr1 pr2 relrisksetx meansetx us_indally minsimqi, prval(3) genpr(pr1)setx us_indally p75simqi, prval(3) genpr(pr2)gen relrisk = pr2/pr1sumqi relriskdrop pr1 pr2 relrisk		 setx meansetx us_indexreltradep1995 0 simqi, prval(2) genpr(start)setx us_indexreltradep1995 75simqi, prval(2) genpr(end)gen diff=end-startsumqi diff, level(95)drop start end diff		 setx meansetx us_indexreltradep1995 0 simqi, prval(3) genpr(start)setx us_indexreltradep1995 75simqi, prval(3) genpr(end)gen diff=end-startsumqi diff, level(95)drop start end diffsetx meansetx us_indexreltradep1995 minsimqi, prval(2) genpr(pr1)setx us_indexreltradep1995 p75simqi, prval(2) genpr(pr2)gen relrisk = pr2/pr1sumqi relriskdrop pr1 pr2 relrisksetx meansetx us_indexreltradep1995 minsimqi, prval(3) genpr(pr1)setx us_indexreltradep1995 p75simqi, prval(3) genpr(pr2)gen relrisk = pr2/pr1sumqi relriskdrop pr1 pr2 relriskdrop b*				*estout a*,  cells(b(star fmt(3)) se(par))   starlevels(* 0.10 ** 0.05 *** 0.01)   stardetach   stats(N, fmt(0) labels("Number of Observations"))   varwidth(18)   modelwidth(8)   mlabels("Model 1" "Model 2" "Model 3" "Model 4" "Model 5" "Model 6")  order(us_SH1995 us_indexreltradep1995 duration geo_salience number_actors coldwar)   varlabels(us_SH1995 "SH Score" duration "Duration" geo_salience "Salience" number_actors "Number of Actors" violence "Violence" us_indexreltradep1995 "Rel. Trade Dep." coldwar "Coldwar"_cons "Constant")  legend   prehead(@title "(Multinomial Logit Estimation, 2 = Diplomatic or Economic Activity, 3 = Indirect of Direct Military Activity)" "" @hline)  posthead(@hline)  prefoot("") postfoot(@hline)  style(fixed)  title("Security and Economic Hierarchy and the Probability of Joining a Crisis, 1950-2000") *estout b*,  cells(b(star fmt(3)) se(par))   starlevels(* 0.10 ** 0.05 *** 0.01)   stardetach   stats(N, fmt(0) labels("Number of Observations"))   varwidth(18)   modelwidth(8)   mlabels("Model 1" "Model 2" "Model 3" "Model 4" "Model 5" "Model 6")  order(us_SH1995 us_EH1995 duration geo_salience number_actors coldwar)   varlabels(us_SH1995 "SH Score" duration "Duration" geo_salience "Salience" number_actors "Number of Actors" violence "Violence" us_EH1995 "EH Score" coldwar "Coldwar"_cons "Constant")  legend   prehead(@title "(Multinomial Logit Estimation, 2 = Diplomatic or Economic Activity, 3 = Indirect of Direct Military Activity)" "" @hline)  posthead(@hline)  prefoot("") postfoot(@hline)  style(fixed)  title("Security and Economic Hierarchy and the Probability of Joining a Crisis, 1950-2000") 